The Stokes Phenomenon and Schwinger Vacuum Pair Production in Time-Dependent 

Laser Pulses 
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Particle production due to external fields (electric, chromo-electric or gravitational) requires evolv- 
ing an initial state through an interaction with a time-dependent background, with the rate being 
computed from a Bogoliubov transformation between the in and out vacua. When the background 
fields have temporal profiles with sub-structure, a semiclassical analysis of this problem confronts 
the full subtlety of the Stokes phenomenon: WKB solutions are only local, while the production rate 
requires global information. Incorporating the Stokes phenomenon, we give a simple quantitative 
explanation of the recently computed [Phys. Rev. Lett. 102, 150404 (2009)] oscillatory momentum 
spectrum of e~ pairs produced from vacuum subjected to a time-dependent electric field with 
sub-cycle laser pulse structure. This approach also explains naturally why for spinor and scalar 
QED these oscillations are out of phase. 

PACS numbers: 12.20.Ds, ll.15.Tk, 03.65.Sq 



The Schwinger effect, the non-perturbative production 
of electron-positron pairs from vacuum in an external 
electric field, is a highly non-trivial prediction of QED 
Q , but the physical scales arc such that it is so weak that 
it has not yet been directly observed. Recent experimen- 
tal advances Q have raised hopes that lasers may achieve 
fields just one or two orders of magnitude below the crit- 
ical field strength of Ecr ~ 10^^ V/cm, either in optical 
high-intensity laser facilities such as HiPER (Rutherford 
Lab.) and the Extreme Light Infrastructure (ELI), or in 
X-ray free electron laser facilities. Theoretically, recent 
analyses suggest that the non-perturbative Schwinger ef- 
fect may be observable at these lower field strengths, by 
careful shaping and combining of laser pulses leading to 
a " dynamically assisted Schwinger mechanism" enhance- 
ment (il-fiol . The most important message is that the de- 
tailed shape of the laser pulse is significant, which moti- 
vates the extension presented here of the standard WKB 
approach to include more realistic laser field profiles. 



Observation of the Schwinger effect in the non- 
perturbative domain has the potential to yield valuable 
insight into analogous gravitational effects UJ, , such 
as Unruh and Hawking radiation, where direct experi- 
ments are not feasible, and where issues such as back- 
reaction and out-of-equilibrium physics are poorly un- 
derstood. The basic physics is also relevant for atomic, 
molecular, astro- and plasma physics with ultra-high in- 
tensity lasers, where non-perturbative effects are crucial 
I3I IJ], for heavy ion collisions 15 1, and for the Landau- 



Zener effect. 

We model the electric field in the focal region of two 
counter-propagating laser pulses by a spatially homoge- 
neous electric field E{t) ~ (0,0, E{t)), with vector po- 
tential A(t) = (0,0, A(i)), such that E{t) = -A{t). Even 
in this approximation where we neglect spatial focussing, 
the laser field may involve many physical scales, leading 
to interesting new phenomena jl6l |. For example, con- 



sider an electric field: 



E{t) = Ea cos{ujt + (f>) cxp {-f/{2T^)) 



(1) 



Here lo is the laser frequency, r defines the pulse length, 
and (j) is the "carrier phase" offset. The first surpris- 
ing result in 16| was that the longitudinal momentum 
spectrum of the produced electron-positron pairs is ex- 
tremely sensitive to the value of wr, even when (/) = 0. 
For LUT > 4, the momentum spectrum exhibits oscilla- 
tions, and these become dramatically enhanced as 4> in- 
creases, to the point where at (f> = 7r/2 the spectrum 
develops minima with zero produced pairs (see Fig. 4 in 
[lij ) . The second surprising result in [3l was that the os- 
cillatory minima and maxima are interchanged between 
spinor and scalar QED. By contrast, when computing 
the total pair production rate (obtained by an integral 
over the momenta), one conventionally approximates the 
case of real spinor QED by scalar QED, with an overall 
multiplicative spin factor of 2. A direct application of 
the usual semiclassical "imaginary time method" (ITM) 



171-12 1| to this problem does not account for these os- 
cillations, let alone for the difference of phase between 
spinor and scalar QED. Here we show that the Stokes 
phenomenon gives a quantitative semiclassical explana- 
tion of both these effects. 

The essential physical interpretation of these oscilla- 
tions is a resonance effect in the corresponding quantum 

This same physical 



16| 



mechanical scattering problem 
explanation has also been noted for the photoclcctron 
spectrum in atomic ionization_ll9i. where such oscilla- 
tions have been observed 



22, 2 



We turn this physical 
picture into a quantitative method. This should also be 



relevant for the matterless double-slit experiment [I6, 24 
Recall that with a time-dependent electric field, the pair 
production process can be reduced to a one-dimensional 
over-the-barrier "quantum mechanical" scattering prob- 
lem 0, [3] , with effective " Schrodinger equation" [in 



2 



t rather than x] 
4> + Q\t)c^^0 



Q^t)^m'+pl + {p-A{t)f (2) 



coming from the Klein-Gordon equation for the particle 
in the presence of the laser field. (We first discuss scalar 
QED and later come to spinor QED, where the relevant 
equation is the Dirac equation.) It is an over-the-barrier 
scattering problem since the "potential", — (p— 
is negative, while the "energy" (m^ + p']_) is positive. 
Thus the reflection coefiicient, from which we deduce the 
probability of pair production, is exponentially small in 
the scmiclassical regime where Eq <C m^. This scat- 
tering problem can be solved using WKB methods, or 
numerically using the quantum kinetic approac h Il6ll25 | 
or direct integration of the scattering problem jlSj . The 
equivalence between these approaches is explai ned in [2 

The Bogoliubov transformation approach 
is based on the field decomposition 




I'Q ^e*^*Q 

2Q" V2Q 



-iQ ( —^=e 



(3) 



which enforces equations relating the coefficient functions 
a and /3. Unitarity requires |ap — = 1, and the parti- 
cle number momentum spectrum is N(j)) = \l3p{t = oo)^, 
related to the reflection coefficient |i?p = |/3p/(l-h|/3p). 
The Stokes phenomenon is relevant because in calculat- 
ing we compare /3{t = oo) to /3{t = — oo). However, 
the leading WKB solutions, e^^ -f / y/2Q , on which ^ 
is based, are multi- valued functions, only deflned locally. 
Evolving a scmiclassical approximation from i = — oo to 
t = +O0, we cross Stokes and anti-Stokes lines, lines along 
which e^* ^ ^ are exponential or oscillatory, respectively. 
On crossing such lines, wc must take care to keep track 
properly of the dominant and sub-dominant solutions. 
This is the Stokes phenomenon 27-30j. 

Suppose the zeros of Q{t), the "turning points" (t.p.'s), 
are first-order. Since this is an over-the-barrier scattering 
problem, the t.p.'s lie off the real axis, in the complex 
plane, and for real laser pulses they occur in complex 
conjugate pairs; furthermore, those closest to the real 
axis tend to dominate in the scmiclassical regime. For 
a simple single-pulse field like E(t) = Escch^ (ujt), with 
A{t) = -E/uta.nh{ujt), or E{t) ^ £'e"("*)', with = 
— ^E I {2uj) Erf(cj<), a single pair of complex conjugate 
t.p.'s dominates, and the phase integral method leads to 
the familiar formula [H, [gl, [1^ [H [si 1 



N{p) 



-2K 



K 



Qdt 



tp 



(4) 



where the integral is along the line joining the two com- 
plex conjugate t.p.'s. In the ITM, one expands in mo- 
menta to obtain the general Gaussian expression [T§| 



N{p) w exp [-S'ci - cipj 



C2P 



(5) 



where Sc\ is the classical action evaluated on the con- 
tour, ci ^-S^S.u and C2 = -2m2g^S'ei. This ITM 



result ([U [18|, [I9j gives a compact expression that is the 
basis for most studies of vacuum pair production, and 
when integrated over momentum to give the total rate 
it gives excellent agreement with numerical (or exact) 
results d, 



18, m 



351 . ([5]) was applied 
However, the ex- 



|21|, |31|, 133. In 

to the envelope field ^ with = 
pression ([5]) clearly cannot exhibit any of the numerically 
observed oscillations, as noted in [16[. This deficiency is 
not cured by higher order terms in the momentum ex- 
pansion in ([5]), nor is it cured by including higher-order 
WKB terms. The problem is that (|4]) and ([5]) are based 
on the assumption of just one pair of turning points, on 
the imaginary axis. However, for complicated fields, with 
sub-cycle structure, the essential shape of the "scatter- 
ing potential", —{p — A{t))'^ , changes dramatically as p 
varies, as illustrated in Fig 1. This can lead to scattering 
resonances, encoded semiclassically in multiple pairs of 
complex conjugate t.p.'s, whose locations are correspond- 
ingly sensitive to variation of p. We show below that the 
oscillatory momentum behavior can be identified with 
interference effects between such pairs of turning points. 
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FIG. 1: The effective scattering potential, V{t) — ~{p — 
A{t))^, for scalar QED, with the field in ©, for three dif- 
ferent values of the longitudinal momentum, in units of m. 
Note that the form of the potential is highly sensitive to the 
value of p. The dashed line denotes the mass level m^. 



To illustrate the oscillatory phenomenon most clearly 
we consider a simple analytic profile field that exhibits 
the effect found in [l^. The maximal oscillation effect 
occurs with carrier-phase </> = 7r/2, so that E{t) is an odd 
function, and A{t) an even function, so we choose 



A{t) 



E/lo 



(6) 



(l-t-a;2 ^2) 

The scattering potential is plotted in Fig. [TJ The alge- 
braic form of the vector potential makes it easy to find the 
turning points, occurring as two complex conjugate pairs: 

t2{pr 



hip) 



Mp)* -i4(p), as 



shown in Fig. 2. Notice that all four t.p.'s are equidis- 
tant from the real axis, for all p. Fig. [5] also shows the 
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Stokes and anti-Stokes lines. Since there are two pairs of 
t.p.'s, the WKB analysis leading to (|4]) and ([5]) must be 
generalized to account for the crossing of multiple Stokes 
and anti-Stokes lines for multiple pairs of t.p.'s in evolv- 
ing from t = — oo to t = +00. This corresponds to the 
case of over-the-barrier scattering with two bumps in the 
scattering potential, which has been solved in |36j using 
the phase integral approximation (PIA). Adapting their 
result, we find the simple expression 
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FIG. 2: The 4 complex turning points, ti, . . . ,t4, for the 
field showing also the anti-Stokes lines (solid, blue, lines). 
Stokes lines (dotted, black, lines) and the integration contours 
(dashed, red, lines) used in ^ and ((8]). 



scalar 



-2Ki I -2K2 



2 cos(2q;) ( 



-Kt-K2 



(7) 





rt2 






Ki = 


/ Qdt 




/ Qdt 




Jti 




Jt3 



L 



a{K) 



L - a{K,) - a{K2) 
ne\ I Qdt 



In fact, in this case Ki = K2, so we can write A'scaiar ~ 
4cos^(a)e~^^i. Note the appearance in ([7]) of the in- 
terference term, cos(2a), involving an integral between 
different pairs of turning points. This term is responsible 
for the oscillations in the momentum spectrum, as shown 
in Fig. [3] where we compare (O with the exact numerical 
result, and with a naive application of the ITM result (j4]), 
just taking the first two terms in ([T]). The agreement of 
([7]) with the numerical result is excellent. A generaliza- 
tion to more than two pairs of t.p.'s is discussed in [STj . 

Having given a quantitative semiclassical explanation 
of the longitudinal momentum oscillations for scalar 
QED, we now turn to spinor QED, for which there is 
a similar scattering formulation 18, [2^, 26 1. The key 



2.x 10 



1.5x10 




1.x 10"'* 



5.x 10 



FIG. 3: Longitudinal momentum spectrum of e pairs for 
the field (O, for £ = 0.2, oj = 0.1, all in units of m. The solid 
lines are our WKB expressions in ((7|) and ((8|, the dashed lines 
are exact numerical results, and the dotted (red) line is the 
naive ITM expression, neglecting the interference term. The 
oscillatory blue lines are scalar QED and the oscillatory black 
lines are spinor QED. The quantitative agreement of ((7| and 
((H]) with the numerics is excellent. 



Bogoliubov coefficients have a reversed sign relative to 
the scalar case: now jap -I- |/3p = 1. This changes the 
form of the F-matrix of the PIA in [3I] , and is ultimately 
related to the double-valuedness of the spinor wavefunc- 
tion. We find the scalar result ([7]) is modified to 



spinor 



-,-2^2 



2 cos(2q!) e 



-Ki-K-2 



(8) 



difference is that the unitarity conditions on the spinor 



where the Ki and a are defined as in ([7|) . The only change 
is the sign of the interference term. (When Ki = ^^2, as 
for the field in wc have TVspinor ~ 4,s\r?{a) e'"^^^ .) 
Physically, this term is produced by interference between 
waves reflected by the double-bump structure, and for 
fermions there is an additional phase shift on reflection, 
which ultimately leads to this sign change. In Fig. |3]we 
plot this spinor result ([8]) and we see that it is in excellent 
agreement with the exact numerical results. The results 
d?]) and dH]) explain clearly why the oscillations are out 
of phase between spinor and scalar QED, and why the 
envelope of the two is the naive ITM result. Of course, 
if one is interested only in the total pair production rate, 
obtained by integrating over p, then the difference be- 
tween spinor and scalar QED is washed out, and agrees 
with the answer obtained by integrating over the enve- 
lope result coming from just the first two terms in ([7]) 
or ([8]), since they oscillate about the same envelope. To 
conclude, we sketch in Fig. 2] the turning points of the 
carrier-phase field in ([T|), for the case = tt/2. Note 
that there are infinitely many pairs of complex conju- 
gate t.p.'s. But the two pairs of t.p.'s closest to the real 
axis dominate, and using formulas (O and ([8]), we re- 
produce the oscillatory behavior of the electron-positron 
longitudinal momentum spectrum found numerically in 
When = the sensitivity to the value of lut can 
also be understood in terms of the location of the t.p.'s. 
Our result is general and simple to use, and has appli- 
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FIG. 4: Contour plot of + (p- for the electric field 
([T]), with = 7r/2, p = and ojr = 4, showing the infinite 
set of pairs of complex turning points. The two central pairs 
closest to the real axis dominate the semiclassical analysis. 



cations beyond this particle-production context, for ex- 
ample to strong-field ionization of atoms and molecules 



13, 19|, to particle production problems in cosmo logy 
I2L 1381 and to quasinormal modes of black holes j39| . 
The basic message is that when the time dependence of 
the external background field has more substructure than 
just a single bump, the usual textbook ITM result (|4]) 
generalizes in an interesting way that requires fuller con- 
sideration of the Stokes line structure of the associated 
scattering problem, and there are important differences 
between the momentum spectra of produced spinor or 
scalar particles. This also has important implications for 
the worldline approach to pair production [40| , which has 
the potential to describe pair production in electric fields 
with both temporal and spatial inhomogeneity. 

We acknowledge the DOE grant DE-FG02-92ER40716, 
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